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*^ ■ We present and solve a model for the susceptibility of two CuC>2 planes cou- 

pled by an interplane coupling Jj_ and use the results to analyze a recent 
"cross-relaxation" NMR experiment on Y2Ba4Cu70is. We deduce that in 
this material the product of J±_ and the maximum value of the in-plane sus- 
ceptibility Xmax varies from approximately 0.2 at T = 200 K to 0.4 at T = 
120 K and that this implies the existence of a temperature dependent in-plane 
spin correlation length. Using estimates of Xmax from the literature we find 
5 meV < J± < 20 meV. We discuss the relation of the NMR results to neutron 
scattering results which have been claimed to imply that in YBa2Cu306+ x the 
two planes of a bilayer are perfectly anticorrelated. We also propose that the 
recently observed 41 meV excitation in YBa2Cu3C>7 is an exciton pulled down 
below the superconducting gap by J±. 
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In the yttrium-barium (Y-Ba) family of high temperature superconductors the basic 
structural unit is a "bilayer", consisting of two CuC>2 planes; the bilayers are separated by 
CuO chains. Neutron scattering ]I[ and more recently NMR experiments || have shown 
that the Cu spins on adjacent planes in a bilayer are coupled. Intra-bilayer coupling has 
been shown theoretically to lead to a "spin gap" |3|-|6| similar to that observed |7| in NMR 
experiments on underdoped members of the YBa family. In view of the great importance 
of the spin gap phenomenon, a quantitative analysis of the spin dynamics of a bilayer is 
desirable. In this letter we provide this analysis and use the results to interpret NMR and 
neutron scattering experiments. 

Focus on the two planes of a bilayer, and neglect coupling to other bilayers. We label 
the spin degrees of freedom by an index a = 1, 2 distinguishing planes, and a site index i. 
We then define the susceptibility 

X a \q^) = | o °°dte^( i? -^)([5f(t),5}(0)]). (1) 

Because of the symmetry under exchange of planes, x ab has only two independent compo- 
nents, x n (g, c<j) = x 22 (<?, ^) and X 12 (g, oj) = x 21 (<L The two independent components 
of x ma y be taken to be the even and odd (under the interchange of planes) components 
Xeven,odd = Xn ±Xi2- In a system with antiferromagnetic coupling one expects that at large 
wavevectors the odd-parity spin fluctuations are softer than the even parity spin fluctuations. 

Further analysis requires a model. We shall assume that the interplane coupling J± is 
weak in comparison to the energies determining the spin susceptibility Xoi.1i w ) °f a single 
plane and that its effects may be modeled via the RPA. Thus we write 

We further assume Xo{Q.i u ) nas the scaling form Xo{Q.i u ) — Xo£ z f where / is a 

scaling function normalized so that /(0, 0) = 1, z is the dynamical exponent, £ is a correlation 
length and q is measured from an ordering wavevector Q which for the present discussion is 
arbitrary. Q is believed to be of the order of (n, it) in high T c materials. From Eq. |] we see 



x- 1 



(2) 



that X U (<1^) = xo(q,u)/(l-(J ±X o(q,uj)) 2 ) and x 12 (g,^>) = J ± Xo(q,u) 2 /(l-(J ± Xo(q,u)) 2 ). 
The crucial parameter controlling the susceptibilities in the static limit is 



We must assume A < 1 so that the material has no long range order. If A 2 1 then 
X U ~ Xo an d x 12 ~ J±Xo- m this limit the interplane coupling has a weak effect and the 
RPA is an appropriate model. On the other hand, if A 2 w 1 then the interplane coupling is 
strong and the use of the RPA may be questioned. 

We now turn to the NMR experiments of interest. These are T 2 experiments performed 
on Y2Ba4Cu70i5, a material in which the single-chain structure of YBa2Cu307 alternates 
with the double-chain structure of YB^C^Os ||. As a result, atoms on different planes of 
a bilayer have somewhat different local environments and therefore somewhat different NMR 
resonance frequencies, which may be independently studied. Despite the differences in local 
environment the electronic properties of the two planes are not very different || , so it is still 
appropriate to model the eletronic properties with Eqs. (1-5). Now the NMR T 2 measures 
the rate at which a nuclear spin is depolarized by interacting with other nuclear spins, i.e. 
it measures the nuclear-spin-nuclear-spin interaction strength. In high T c materials the 
dominant contribution to the nuclear-spin nuclear-spin interaction comes from polarization 
of electronic spins, and may be related to the static limit of the real part of the electronic 
spin susceptibility. In Y 2 Ba4Cu70i5 it is possible to measure T 2 , the rate at which a spin 
in one plane is depolarized by spins in the same plane, and T 2 j_, the rate at which a spin 
in one plane is depolarized by spins in the other plane. T 2 is related to the electronic spin 
susceptibility by || 



A = J ±X 



max 



J±X^ 2 - 



(3) 





while T 2 j_ is given by [Kj 
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We have calculated T 2 and T 2 ± from Eqs. (^],f|,|5]). The precise values obtained depend 
upon the form chosen for f(q£). We have used two forms for f(x) = f(q£,u) = 0): a 
Lorentzian, f(x) = 1/(1 + x 2 ), and a Gaussian, f(x) = exp(— log(2)a; 2 ) ( the log(2) is 
introduced so f(x = 1) = 1/2). We measure £ in units of the lattice constant, we set h — 1 
and assume that the hyperfine coupling can be approximated by its value at Q, Aq. We 
find 



^ = ^xo^m(A) (6) 



T, 
and 

1 



A 2 Q J ±x teg±(A) (7) 



T, 



2± 



where gi n and g± are defined in terms of the function f(x) via 



<4(A) = -j- 
Ztt 



(ia; . „ „„ NO — \ dx 



— (i _ A 2 / 2 ) 2 2vre 2 1^0 1 - A 2 / 2 



(8) 
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In writing Eqs. (8,9) we have assumed that the correlation length is so long that lattice 
effects may be neglected. We have investigated this issue by performing the exact integrals 
numerically. The parameter governing the size of the lattice effects is (7r£) -1 ; for £ > 1 we 
have found that they are negligible. 

In Fig. 1 we present the calculated results for T 2 /T 2 j_ = Agj_(A)/gj n (A). From the 
experimental values T 2 /T 2± = 0.15 at T = 200 K and T 2 /T 2± = 0.30 at T = 120 K g] we 
obtain A w 0.2 at T = 200 K and A w 0.4 at T = 120K. Note that even at the lowest 
temperature we find that the system is in the small A regime in which T 2 /T 2 ± is linear in 
A, suggesting that the between-planes coupling is sufficiently small that the RPA formula 
is justified. 

Now Y 2 Ba 4 Cu70i 5 has a doping somewhere between YBa 2 Cu 3 07 and YBa 2 Cu 3 6 .7 
It seems clear that the strength of the magnetic correlations is a relatively rapid function 



of doping and increases as one moves from O-j to the insulator, while J± is unlikely to be a 
sensitive function of doping. It therefore seems very likely that for dopings corresponding to 
YBa 2 Cu306.7 and below, Xmax is so large that A s=s 1 and the planes are so strongly coupled 
that the properties are not linear in J±. 

The magnitude of J± may be determined if Xmax is known and conversely. If the suscep- 
tibility were only weakly g-dependent then the measured uniform susceptibility Xuniform ~ 
2 states/eV-Cu [|ll|] would provide a good estimate for Xmax and our value A rj 0.4 would 
imply J± ps 0.2 eV. Such a value is very difficult to justify on microscopic grounds because 
the insulating antiferromagnetic parent compounds of the high T c superconductors have in- 
plane exchange constants Ji n - P iane ~ 0.12 eV and it is generally believed that J± <C J% n ~piane- 
Therefore, we believe the cross- relaxation results imply Xmax 3> Xuniform- A similar con- 
clusion has been drawn from an analysis of the magnitude of the in-plane T 2 measured on 
YBa 2 Cu307 0, combined with various assumptions about magnitudes of hyperfine cou- 
plings. The magnitude of the hyperfine coupling drops out of the present analysis. The 
estimate J± ~ 10 — 20 meV has been obtained from band structure calculations |I2[], im- 



plying Xmax ~ 40 states/eV-Cu. The in-plane T 2 experiment led to the value Xmax ~ 80 
states/eV-Cu implying Jj_ ps 5 meV. 

In summary, the cross-relaxation experiment shows that the real part of the susceptibil- 
ity at some non-zero q is much larger than the uniform susceptibility. Now the temperature 
dependence of the T 2 rates must be due to the temperature dependence of this antiferromag- 
netic maximum. Two scenarios have been proposed for the temperature dependence: in the 
antiferromagnetic scenario the temperature dependent quantity is the correlation length 
£. In the generalized marginal fermi liquid scenario the temperature dependent quantity 
is the overall amplitude x PI - From Eqs. (0,0) we see that in the regime where T 2 j_ 
is linear in A the antiferromagnetic scenario predicts T$/T 2 ± is temperature independent, 
while the marginal fermi liquid scenario predicts T 2 2 /T 2 _|_ is temperature independent. The 
experimentally determined ratios are plotted in Fig. 2 and are more consistent with the 
antiferromagnetic scenario. 



The imaginary parts of the two independent susceptibilities Xeven and Xodd are measur- 
able via neutron scattering because they have different dependences on q z , the momentum 
transverse to the CuC>2 planes |14]]. Neutron scattering experiments have been performed 
on a variety of members of the yttrium-barium family of high-T c materials |Il,|I~4}4T6[j . The 



experimental result is that only Xodd is seen. At frequencies less than 30 meV and tempera- 
tures less than room temperature the even parity fluctuations are claimed to be completely 
frozen out. The theory of neutron scattering in high T c materials is presently controversial. 
There is no generally accepted model which correctly accounts for the observed lineshapes 
and temperature dependences. To investigate the connection between the cross-relaxation 
experiments and neutron scattering we have chosen to calculate the ratio of the g-integrated 
even and odd parity susceptibilities. This ratio is insensitive to the precise details of the 
susceptibilities. For definiteness we used the "MMP", dynamical exponent z = 2 ansatz 
Xo(q,uj) = x/(£~ 2 + q 2 — iuj/T). Here T is a microscopic spin relaxation time. The results 
depend on A and on ujsf — r/£ 2 , which is the softest spin fluctuation frequency of a single 
plane. Of course J± will reduce this frequency for the odd parity channel and increase it 
for the even channel. Results are shown in Fig. 3 for several values of A. We see that the 
relative weight of the even parity fluctuations becomes small only for A > 0.5. We believe 
that the neutron results, which seem to require a A > 0.5, are not in contradiction to our 
analysis of the cross-relaxation experiment, which yielded a A < 0.4, because the strongest 
neutron evidence for locked bilayers was obtained from a study of YB^CuaOg.s 0, which 
as we have previously noted is closer to the magnetic instability than Y 2 Ba 4 Cu70i5, and 
therefore may be expected to have a larger A. 

Another experiment in which the even parity fluctuations were not seen at all was an 
observation of a rather sharp peak at an energy of 41 meV in the superconducting state 
of YB^CusCV In this material there is some evidence of a spin fluctuation peak in the 
normal state at a similar energy, but it is much broader in q and uj. We suggest that the 
41 meV peak is an exciton pulled down below the superconducting gap edge by the interplane 
coupling Jj_. We note that YE^CuaOy has somewhat weaker magnetic correlations than 



Y 2 Ba 4 Cu70i5, so one would expect A < 0.4 in the normal state. Now in the superconducting 
state, the presence of the gap implies that x" becomes very small at energies less than the 
gap, and has a peak (the details of which depend in a complicated way on the details of 
the superconducting order parameter) at the gap edge. This behavior of x" implies via 
the Kramers-Kronig relation a large enhancement of the real part, x\ a t the gap edge. We 
propose that this enhanced x leads to a weakly damped pole, at a frequency of order the gap 
edge, in the odd parity channel of the RPA formula, Eq. ^|, and that this pole produces the 
feature seen experimentally. Our proposal provides a natural explanation for the sharpness 
of the feature and for the nearly perfect bilayer correlation observed at low temperatures. 
If this scenario is correct the even parity component of the neutron scattering signal should 
appear above T c . 

In summary, our analysis of the cross-relaxation experiment suggests that the between- 
planes coupling J± is small, but has non-negligible effects which furthermore increase as the 
temperature is decreased. These are precisely the assumptions made in the theories which 
attribute spin gap formation to interplane pairing, so we believe our results tend to support 
this picture. 

We thank T. M. Rice, R. Stern and M. Mali for helpful discussions. H. M. acknowledges 
the hospitality of AT&T Bell Laboratories where part of this work was done. 
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FIGURES 




FIG. 1. Ratio of cross-relaxation rate l/T2± to in-plane relaxation rate I/T2 plotted versus 
coupling parameter A = J±Xmax for Lorentzian (dotted line) and Gaussian (solid line) form factors 
and calculated from Eqs. (6-9). 
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FIG. 2. Experimentally determined ratio of 1/T2± to nth power of I/T2 for n=l (□), 2 (•), 3 
(o) in arbitrary units. That the n=3 (o) curve has less temperature dependence than the n=2 (•) 
curve suggests the existence of a growing magnetic correlation length. 
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FIG. 3. Calculated ratio of (/-integrated odd-parity neutron absorption to q-integrated even 
parity neutron absorption, plotted versus frequency for A — J iSK,max — 

0.0, 0.1 ...0.9. A = 0.0 

corresponds to the lowest curve and A = 0.9 to the top curve. 
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